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Abstract. We model the lOyum absorption spectra of nonspherical particles composed of amorphous silicate. We consider two 
classes of particles, compact ones and fractal aggregates composed of homogeneous spheres. For the compact particles we con- 
sider Gaussian random spheres with various degrees of non-sphericity. For the fractal aggregates we compute the absorption 
spectra for various fractal dimensions. The 10 pm spectra are computed for ensembles of these particles in random orienta- 
tion using the well-known Discrete Dipole Approximation. We compare our results to spectra obtained when using volume 
equivalent homogeneous spheres and to those computed using a porous sphere approximation. We conclude that, in general, 
nonspherical particles show a spectral signature that is similar to that of homogeneous spheres with a smaller material volume. 
This effect is overestimated when approximating the particles by porous spheres with the same volume filling fraction. For ag- 
gregates with fractal dimensions typically predicted for cosmic dust, we show that the spectral signature characteristic of very 
small homogeneous spheres (with a volume equivalent radius r v < 0.5 fim) can be detected even in very large particles. We 
conclude that particle sizes are underestimated when using homogeneous spheres to model the emission spectra of astronomical 
sources. In contrast, the particle sizes are severely overestimated when using equivalent porous spheres to fit observations of 
10 pm silicate emission. 

Key words. Infrared: general - Stars: circumstellar matter - Stars: planetary systems: protoplanetary disks 



1. Introduction 

The interpretation of absorption and emission spectra observed 
from astronomical objects requires knowledge of the absorp- 
tion cross section as a function of the dust grain characteris- 
tics, such as size, shape and composition. Usually the absorp- 
tion spectra are modeled using homogeneous spherical parti- 
cles for which c alculations can easily be done using Mie the- 
ory jMieL Il908h . Although cosmic dust grains are in general 
not homogeneous spheres, in some cases Mie theory calcu- 
latio ns can reproduce the observations quite accu rately (see 
e.g. lHansen & Hovenieiifl974l:lKemper et all|2004l) . However, 
in other cases one has to find a way of modeling the effects 
of particle shape in order to r eproduce observations or labo- 
ratory measurements (see e.g. Mishchenko et al., 2000). It is 
therefore important to know the effects of the adopted particle 
shape model on the derived dust parameters, such as the par- 
ticle size and structure (see e.g. lFabian et allfeoOltlMin et all 
l2003[l2005al) . 

From different formation mechanisms different types of 
particles may form. For example, when dust grains form from 
direct gas phase condensation, compact particles may be cre- 
ated. Alternatively, when dust grains stick together to form 
larger particles, complex aggregated structures may be formed. 
We study particles in both classes using irregularly shaped 



compact particles and fractal aggregates. For the compact 
^article shapes we use so-called Gaussian random spheres 
( Muinonen SLafl 1996) which allows for a varying degree of 
irregularity. The fractal aggregates are composed of homoge- 
neous spheres and are constructed using a sequential tunable 
particle cluster aggregation method, that allows us to construct 
particles with arbitrary values of the fractal dimension. 

We concentrate on the effects of particle non-sphericity on 
the so-called 10/iin amorphous silicate feature. Amorphous sil- 
icate with an olivine-type stoichiometry is one of the most 
abundant dust components in various astronomical environ- 
ments. We will use this type of amorphous silicate and in the 
rest of this paper refer to it as amorphous olivine-type silicate. 
The absorption spectra of small amorphous silicate grains dis- 
play a characteristic feature peaking around 10 pm. The wave- 
length position of this feature is ideal for ground based obser- 
vations because of the 10 /mi atmospheric window. 

The 10 pm absorption spectra are c omputed using the 
Discre te Dipole Approximation (DDA). Ha ge & Green berg 
(1990) proposed to approximate the optical properties of fluffy 
aggregates with those of so-called equivalent porous spheres. 
In their method they take a porous sphere with the same vol- 
ume filling fraction as the fluffy aggregate and compute the op- 
tical properties of this porous sphere using effective medium 
theory. Using this approximation they concluded that the spec- 
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tral signature of fluffy aggregates is similar to that of compact 
spheres with a much smaller volume. We test the validity of this 
approach for calculations of fractal aggregates and complex 
shaped compact particles by comparing the spectra resulting 
from this approximate method with those obtained using DDA 
calculations. In addition we compare the spectra obtained with 
DDA to those of volume equivalent homogeneous spheres and 
analyze the error that is made when employing homogeneous 
spheres to fit observations of the emission spectra of nonspher- 
ical particles. 

The method is briefly outlined in Sect. [2] The particle 
shapes are described in Sect. [3] In Sect. @] we present the re- 
sulting lOyimi amorphous silicate spectra. An analysis of our 
results and a discussion of the implications on the modeling of 
10 fim spectra in astronomical observations is given in Sect. [5] 

2. Absorption and emission spectra 

The shape of the absorption spectrum of a dust grain, i.e. its 
absorption cross section as a function of wavelength, contains 
important information about the dust characteristics. This spec- 
trum can be observed as absorption against a strong infrared 
background, or as thermal emission. 

In this paper we consider two measures for the size of a 
particle, the volume equivalent radius and the circumscribed 
sphere radius. The volume equivalent radius, ry, is defined as 
the radius of a sphere with the same material volume as the par- 
ticle. For very fluffy particles, i.e. particles with a small space 
filling fraction, the volume equivalent radius does not provide 
an accurate measure for the actual linear extent of the particle. 
For this we define the circumscribed sphere radius, r c , to be the 
radius of the smallest sphere centered on the center of mass of 
the particle and containing the entire particle. As an important 
parameter we consider the ratio of the circumscribed sphere 
radius and the volume equivalent radius 



r = 



(1) 



(2) 



This parameter determines how dense the material is packed in 
space. A very fluffy aggregate will have a high value of y while 
a very dense particle will have a value of y close to unity. For a 
homogeneous sphere y — 1 . 

The mass absorption coefficient, k, of a dust grain is defined 
as its absorption cross section per unit mass 

_ ^ a t>s 

K ~ ~M~ ' 

Here C a bs is the absorption cross section, and M is the mass of 
the particle. For grains that are much smaller than the wave- 
length of radiation both inside and outside the particle, i.e. in 
the Rayleigh domain, the absorption cross section is propor- 
tional to the particle material volume. This means that the mass 
absorption coefficient of very small particles, which is gen- 
erally relatively high, is independent of the particle size. For 
compact particles much larger than the wavelength we are in 
the geometrical optics domain and the absorption cross section 
is proportional to the surface area of the particle. For a compact 
particle this implies that the mass absorption coefficient scales 
as ry l for large values of ry. 



In general the optical properties of nonspherical particles 
depend on the orientation of the particle with respect to the 
incoming electromagnetic field. In the astronomical environ- 
ments we are interested in, we are dealing with an ensemble 
of dust grains with random orientations. Therefore, in this pa- 
per we will always average the absorption cross section over all 
particle orientations. 

2.1. The Discrete Dipole Approximation 

The discrete Dipole Approximation (DDA) proposed by 
IPurcell & Pennvpackerl Jl973t) allows for the computation of 
the optical properties of arbitrarily shaped particles. In the 
DDA the particle material volume may be divided into small 
volume elements. Each of these volume elements is assumed 
to interact with the radiation as a single dipole. The interac- 
tion of all dipoles with the incoming field and with the field of 
all other dipoles is then obtained by solving a 3N x 3N matrix 
equation, where N is the number of dipoles used to represent 
the particle volume. For a theoretical foundation of the me thod 
in terms of the Maxwell equa tions see Lakhtajdl dl992l) and 
lLakhtakia & Mulhollandl il993b . For applications of the DDA 
in the Rayleigh limit, i.e. for sizes of the particles as a w hole 
small c ompared to the waveleng th, see lHenning & Stognienkol 
dl993l) : IStognienko et alJdl995l) . 

In order for the DDA to be valid, the size of the vol- 
ume elements has to be much smaller than the wavelength 
of radiat ion both inside and outside the particle, i.e. (see 
e.g^Draine & Goodmanlll993rlDraine & Fl atau. 1994; Draine, 



ka\m\ < 1, 



(3) 



where k - 2n/A is the wavenumber of radiation in vacuum, A 
the wavelength in vacuum, a the size of the volume element, 
and m the complex refractive index of the material. 

DDA computations for large particles require a large num- 
ber of dipoles. Several methods have been suggested to speed 
up the computations. The most important one makes use 
of a Fast Fourier Transform (FFT) t ogether with a conju- 
gate gradient solution method (see e.g. lGoodman et all 1 199 it 
iHoekstra et all 119981) . A disadvantage of this method is that 
all dipoles have to be located on a rectangular grid, conse- 
quently also the empty, vacuum places on the grid are taken 
into account in the computation. This slows down the method 
for fluffy particles with a low volume filling fraction. For ex- 
ample, for particles with a volume filling fraction of 1%, N 
is increased by a factor 100. Therefore, we chose for a direct 
solution method of the DDA equations. For moderate values 
of N this allows for relatively fast computations, not only for 
compact but also for very fluffy grains. An additional advan- 
tage is that we only have to do the computation for a particle 
in a fixed orientation once in order to obtain the optical prop- 
erties of the particle for any other orientation, whereas the FFT 
method requires the computation to be repeated for every ori- 
entation. Since the computing time required for the direct so- 
lution method scales as A^ 3 , and the method using FFT scales 
as N log N, the FFT method becomes faster when large values 
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of N are needed. Typically, both methods are equally fast for 
N ~ 10 3 for a particle with a volume filling fraction of 1%. 



spheres, while for the aggregates we use fractal aggregates with 
varying fractal dimensions. 



2.2. Porous spheres 

Another method that is frequently employed to obtain 
the absorption and emission spectra of aggregated par- 
ticles is that of approxima tin g the particles b y porous 
spheres faage & Greenbergl Il990t iMin et all l2005al 
IVosh chinniko v et all 120051) . The parameter y can be eas- 
ily transformed to the so-called porosity factor, P, which 
is defined as the vacuum fraction of the volume inside the 
circumscribing sphere and is given by 



1 



= 1 - 



(4) 



lHage & Greenbergl Jl99dh use the porosity factor to describe 
the optical properties of fluffy aggregates. In their approach 
they take a solid sphere with radius r c and an effective re- 
fractive index depending on the porosity factor. In this way 
a porous sphere is created with the same volume equiva- 
lent radius as well as the same circumscribing sphere radius 
as the original particle. We will refer to these particles as 
equivalent porous spheres. To comp ute the effective refrac- 
tive index Hage & Greenberg ( 1990) apply the Garnett mix- 
ing rule (Maxwell-Garnett, 1904, usually refered to as the 
Maxwell-Garnett mixing rule), taking the material as inclu- 
sions in a vacuum matrix. Then the effective complex refrac- 
tive index, , of a partic le with porosity factor P is given by 
( lHage & G reenber gj Il99dl) 



'eff 



1 + 



3(1 -P)(m 2 - l)/(w 2 + 2) 
1 -(1 -P)(m 2 - l)/(m 2 + 2)' 



(5) 



with m the complex refractive index of the material. 

The combination of Garnett effective medium approxima- 
tion with exact Mie calculations for volume equivalent spheres 
will be referred to as the porous sphere approximation through- 
out this paper. The validity of this approximation was tested by 
lHage & Greenbergl Jl990l) for particles where the volume ele- 
ments are much smaller than the wavelength of incident radia- 
tion and in addition are randomly distributed over the circum- 
scribing sphere volume. They found that in this case, the porous 
sphere approximation yields rather accurate results. However, 
in realistic aggregates the constituents are not randomly dis- 
tributed in space. In contrast, the constituents have to be touch- 
ing to form, for example, fractal like structures. We will test the 
validity of the porous sphere approximation for computations 
of emission spectra of fractal aggregates and complex shaped 
compact particles by comparing the spectra resulting from this 
approximation with those obtained from DDA computations. 

3. Particle shapes 

The shape of a dust grain is an important parameter determin- 
ing its spectroscopic characteristics. In this paper we distin- 
guish between two classes of shapes, aggregates and compact 
particles. For the compact particles we use Gaussian random 



3.1. Gaussian Random Spheres 

A Gaussian random sphere Jlvluinonen et all[l996l) is a homo- 
geneous sphere the surface of which is distorted according to 
a Gaussian random distribution. The distortion of the surface 
is parameterized by two shape variables, the standard deviation 
of the distance to the center, cr, and the average correlation an- 
gle, T. Using various values for these two parameters allows us 
to construct particles with varying degrees of irregularity. The 
number of 'hills and valleys' on the surface in a solid angle is 
determined by the value of T, while cr determine s the height o f 
these hills and valleys. For details see Muin onen et all l fl99rj|) . 

In Fig.[flwe show pictures of the Gaussian random spheres 
used in our computations. We vary the value of cr from 0.1 to 
0.7 and fix the value of T to 10°. 

3.2. Fractal aggregates 

In astronomical environments large grains form by aggrega- 
tion of small particles (see e.g. Kempf et all 1 19991) . Studies on 
interplanetary dust particles from likely cometary origin show 
that cometary dust grains are likely to be aggregates. A special 
class of aggregates are the so-called fractal aggregates. 

A fractal aggregat e composed of hom ogeneous spheres 
obeys the scaling law (Fili ppovetaTll2000l) 



N = k f I — 



(6) 



Here N is the number of constituents each with radius a; kf is 
the fractal prefactor; Df is the fractal dimension, and R g is the 
radius of gyration defined by 



l N 



(7) 
(8) 



where r, is the position of the z'th constituent, and r () is the cen- 
ter of mass. The value of the fractal dimension can in theory 
vary between the two extremes Df = 1 (a thin, straight chain 
of particles) and Df = 3 (a homogeneous sphere). 

We use a sequential tunabl e par t icle- cluster aggregation 
method developed by Fili DPOv et all J2000h to construct the 
fractal aggregates. This method allows us to construct aggre- 
gates with arbitrary fractal dimension. 

Note that for a fractal aggregate the circumscribed sphere 
radius and the volume equivalent radius cannot be the same. 
The reason is that the constituents are homogeneous spheres, 
therefore they do not allow for a full packing of the volume 
of the circumscribed sphere. The maximum packing densi ty of 
a col lection of spheres with the same radius is n/ VT8 (Hales, 
119921) . This implies that for an aggregate composed of homoge- 
neous spheres the minimum value of y = ( Vl8/7r) I/3 * 1.105. 
The maximum value of y depends on the number of spheres 
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Fig. 1. Pictures of the Gaussian random spheres employed in our calculations. The diameters of the corresponding volume 
equivalent spheres are indicated by bars below the particles. 



that make up the aggregate. The least dense packing is a straight 
line of spheres and results in y = N 2 ^ . 

In nature fractal aggregates can grow with a wide range of 
fractal dimensions depending on the formation mechanism and 
the environment in which they form. In astronomical environ- 
ments the fractal dimension can vary between approximately 1 
up to almost 3. Almost linear aggregates {Df « 1) might form 
when the particles h ave an electric or magnetic dipole and are 
in an external field (Dominik & Nubold, 2002; Nubold et al., 
2003|. In the case that grain growth occurs by coagulation of 
aggregates with approximately equal size the resulting aggre- 
gates will have a fr actal dimension of approximately 1.8 to 2.1 
iKenrof etallll999fc . On the other hand, when the growth oc- 
curs by adding single monomers to a larger aggreg ate, aggre- 
gates with a fractal dimension of 3 might form (Ball & Witten, 

In this paper we consider aggregates composed of homo- 
geneous spheres with fractal dimensions ranging from 1.2 up 
to 2.8. The radius of the monomers, a, is chosen to be 0.4 /mi. 
For particles composed of an amorphous silicate, the absolute 
value of the refractive index in the 10 /mi region reaches val- 
ues up to \m\ ~ 2.2 with a corresponding maximum value of 
ka\m\ = 0.5. Thus, when applying DDA in this spectral region 
for the material we consider, each monomer can be represented 
by a single dipole (see Eq. [3Jl. By using several values of kt 
we find that the results are not very sensitive to the fractal pref- 
actor chosen. For all aggregates we use kf = 2.0. Pictures of 
the fractal aggregates we constructed for various numbers of 
monomers, and thus for various volume equivalent radii, are 
shown in Fig.|2] The bar below the pictures shows the diame- 
ter of a volume equivalent sphere. Also denoted are the corre- 
sponding values of y. Thus we find that for a fractal aggregate 
with Df = 1.2 and a volume equivalent radius ry = 6 /mi, we 
have r c = y ry ~ 320 /mi. 

When we consider growth, we do not scale the aggregate 
but add more grains to it. This is an important difference with 
the Gaussian random spheres, since it implies that aggregates 
with different sizes have different shapes. For Gaussian random 
spheres the size is changed by scaling the entire grain, which 
implies that the shape of the particle is independent of its size. 



4. Results 

4.1. The 10 fim amorphous silicate feature 

The absorption spectrum of silicate particles display a feature 
around 10 /mi due to the stretching mode of the Si-0 bond in 
the Si04 tetrahedra of which silicates are made. Here we will 
concentrate on the 10 /mi feature of amorphous olivine-type sil- 
icate, which is the most important dust constituent in many as- 
tronomical environments. The chemical formula of amorphous 
olivine-type silicate is Mg2. v Fe2-2. v Si04, where < x > 1 gives 
the iron to magnesium ratio. We use t he refractive index as a 
function of wavelength as measured bv lPorschner et all {f995 ), 
for olivine with x * 0.5. 

The spectral shape and amplitude of the mass absorption 
coefficient of particles composed of an amorphous silicate as 
a function of wavelength in the 10 fim region is sensitive to 
the size and shape of the grains causing it. Computations us- 
ing homogeneous spheres show that an increase in particle size 
tends to broaden and flatten the 10 fim feature. When Mie the- 
ory is employed to compute the emission spectrum, the feature 
strength diminishes rapidly when the grai n size is increased. 
Using the porous sphere approximation Hage & Greenberg 
( 1990) have shown that this change in feature strength and 
shape is strongly affected by the porosity factor of the particle. 
Using approximate methods to compute the emission spectra 
they showed that the spectra of large, porous spheres resemble 
those of small compact ones. 

4.2. Fractal aggregates versus Gaussian random 
spheres 

The solid curves in Fig.[3]display the mass absorption coeffi- 
cients as a function of wavelength for the various Gaussian ran- 
dom spheres (Fig.^, and different particle sizes. Also shown in 
this figure is the mass absorption coefficient of volume equiv- 
alent homogeneous spherical particles (dotted lines), and the 
mass absorption coefficient of equivalent porous spheres with 
porosity factor given by Eq. @ (dashed lines). It is clear that 
in all cases an increase of the particle size beyond 1 /mi, leads 
to a decrease of the strength of the 10 fim feature. Also, for the 
grains with a volume equivalent radius larger than 1 /mi, in- 
creasing the value of cr increases the peak value. The reason 
is probably that the bumps and spikes sticking out of the parti- 
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Fig. 2. Pictures of the fractal aggregates we used in our calculations. Indicated by the bar below the aggregates is the diameter of 
the corresponding volume equivalent sphere. 



cle for high values of cr (see Fig.^ roughly interact with the 
radiation like small separate particles. In general, the feature 
as computed using volume equivalent homogeneous spheres is 
narrower and weaker. Also, it is clear that the feature computed 
using the porous sphere approximation is in general narrower 
and stronger, thus resembling the feature of much smaller par- 
ticles. When for homogeneous spheres the grain size is in- 
creased, the feature weakens and flattens, i.e. it becomes more 
flat-topped. For the Gaussian random spheres the feature also 



weakens when the grain size is increased, however, the shape 
of the feature still shows a more rounded top. 

By considering the spectra of different realizations of the 
Gaussian random spheres with the same shape parameters cr 
and r but different seeds for the random number generator we 
find that the differences in the absorption spectra are small and 
that all trends are conserved (not shown in the figure). 

Plotted by the solid curves in Fig. |4] is the mass absorp- 
tion coefficient as a function of wavelength for the fractal ag- 
gregates with different fractal dimensions and aggregate sizes 
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Fig. 3. The mass absorption coefficient, k, averaged over all particle orientations as a function of wavelength for Gaussian Random 
Spheres composed of amorphous olivine-type silicate (solid lines) with various values of the shape parameter cr and grain sizes 
ry. Also shown are the spectra for volume equivalent homogeneous spheres (dotted lines) and for equivalent porous spheres with 
a porosity factor given by Eq. (|4j (dashed lines). Particle size is increased from top to bottom, while the particle non-sphericity is 
increased from left to right. Note that in accordance with the shape of the particles, the porosity of the grains increases from left 
to rig ht. The refractive index of amorphous olivine-type silicate as a function of wavelength used in the calculations was taken 
from lDorschner et all (1 19951) . 



as shown in Fig. |2] It is clear from this figure that the 10/mi 
feature of fractal aggregates with a low value of the fractal 
dimension in general has a much weaker dependence on the 
particle size than that of aggregates with a higher value of the 
fractal dimension. For a fractal aggregate with Dj = 1.2 the 
feature strength is only very slightly decreased even if the vol- 
ume equivalent radius is increased from ry = 0.4 /vm to 6^m, 
in which case r c = 320pm. The reason is probably that the 
constituents in fractal aggregates with low fractal dimensions 
are separated by large distances. This reduces the interaction 
between the different volume elements within the aggregate. 
Thus, the monomers interact, to a certain extent, as if they were 
separate small particles, displaying a feature typical for small 
dust grains. We have to note here that there is a difference be- 
tween the spectra for the separate constituent (upper panels) 



and the fully grown aggregates (shown in the lower panels) 
even for the most fluffy aggregates. This means that even in the 
most extreme case Df = 1 .2 there has to be interaction between 
the constituents of the aggregate. Thus, the effects of aggrega- 
tion on the spectral signature will always be visible, even in the 
case that the aggregates form almost linear chains. In general, 
when increasing the fractal dimension, the spectral signature 
of the monomers becomes more apparent in the spectrum of 
the aggregate. Comparing the resulting spectra for the fractal 
aggregates with homogeneous spheres we see that, as with the 
Gaussian random spheres, the features for the larger fractal ag- 
gregates also have a generally more rounded shape. The spectra 
of very compact aggregates (Df = 2.8) show virtually no dif- 
ferences with those of homogeneous spheres. 
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Fig. 4. Same as Fig.[3]but for fractal aggregates with various values of the fractal dimension (solid lines). The increase in particle 
size for aggregates implies an increase of the number of constituents (see Fig.|2ji. Note that in this figure the porosity increases 
not only from left to right but also from top to bottom (see also Fig.|3 



Also for the fractal aggregates we considered different re- 
alizations of the particles with the same shape parameters but 
different seeds of the random number generator and found that 
the differences in the spectra obtained are negligible. 

Comparing the features obtained using DDA calculations 
for fractal aggregates with those obtained using the equiva- 
lent porous spheres (dashed curves in Fig. |4j it is clear that 
the effect of fluffiness is overestimated using porous spheres. 
The features obtained using this approximation are, in general, 
much sharper, stronger and narrower than those obtained using 
DDA calculations, thus resembling more the feature of the con- 
stituents of the aggregates. The effective medium theory used 
to compute the spectra of porous spheres assumes that the con- 
stituent particles of the aggregate are distributed randomly over 
the entire volume of the circumscribing sphere. As mentioned 
before, this is in sharp contrast with realistic aggregates where 
the constituents have to be in contact, thus increasing the in- 
teractions between the aggregate constituents. In this way it is 
easy to understand that the porous sphere approximation un- 



derestimates the effects of grain growth for fluffy aggregates. 
Our detailed calculatio ns qua litatively confirm the conclusion 
by Hage & Greenberg ( 1990) that by increasing particle fluffi- 
ness the spectral signature of the constituent particles of the 
aggregate become more visible. However, we also have to con- 
clude that the porous sphere approximation severely overesti- 
mates this effect for realistic aggregates. 

5. Implications for grain growth and dust modeling 

In this section we will quantify the effects of the assumption 
of particle shape on the derived grain size and dust mass when 
fitting observed infrared lOyum spectra. We do this by taking 
the lOyum spectra presented in the previous section and fitting 
these with the 10 fim spectrum of an ensemble of homogeneous 
spheres with radius rj and total mass Mf. In the same way we 
fitted the spectra using an ensemble of porous spheres, apply- 
ing the porous sphere approximation, with volume equivalent 
radius r p and total mass M p . The porosity of these spheres was 



8 



M. Min et al.: The 10//m amorphous silicate feature of fractal aggregates and compact particles with complex shapes 



taken to be equal to the porosity of the particle whose spec- 
trum was fitted. The fits using homogeneous spheres showed 
that, in general, the spectral shapes shown in Figs. 0] and [3] are 
very well reproduced using homogeneous spheres. In almost 
all cases the fitted curves deviate less than the linewidth from 
the curves shown in Figs. [3] and 0] Only for the large (4 and 
6 //m) Gaussian random spheres with high values of cr and the 
large fractal aggregates with intermediate values of Df the fits 
are of a slightly lesser quality. In these cases, the nonspherical 
particles display a feature with low contrast combined with a 
relatively rounded top. This cannot be reproduced accurately 
by using homogeneous spherical particles. The fits using the 
porous sphere approximation showed also quite good agree- 
ment. However, for the large (r v > 2 fim), very fluffy aggre- 
gates (Df = 1.2 and 1.8), the agreement is much less satisfy- 
ing. The fact that most spectra are fitted very well implies that 
it is difficult to determine particles shapes from the shape of 
the 10/mi feature alone. Thus, one has to resort to other ob- 
servables, like for example the degree of linear polarization in 
the visible part of the spectrum, to obta in information on the 
particle shape (see e.g. lMin et alll20Q5hh . 

The resulting values of the best fitting rj and r p for the dif- 
ferent values of cr, Df, and ry are presented in Tableland the 
corresponding values of Mf and M p in Table[2] In Fig.[5]we plot 
rf as a function of ry for all particle shapes we considered. It is 
clear that in almost all cases the size of the particles is underes- 
timated when fitting the 10 fim absorption spectrum using ho- 
mogeneous spheres. This is a consequence of the fact that when 
using homogeneous spheres the particle mass is concentrated in 
a much smaller volume than when realistically shaped particles 
are used. This increases the interactions between the volume 
elements of the particle, leading to a stronger effect of particle 
size on the 10 fim feature when using homogeneous spheres. In 
contrast, when using equivalent porous spheres to fit the spec- 
tra, the particle size is severely overestimated. This in turn is 
caused by the fact that when using equivalent porous spheres, 
the mass is distributed over a much larger volume than when 
using realistic particle shapes, thus leading to a much weaker 
effect of particle size. For the fractal aggregates with very low 
fractal dimension these effects are most extreme. Extrapolation 
of the curves in Fig.|5]beyond 6 fim for fractal aggregates with 
low fractal dimension suggests that even for very large aggre- 
gates the 10 fim feature will display a signature typical for rel- 
atively small homogeneous spheres. The mass estimate we get 
from the fit using homogeneous spheres is in most cases close 
to the real mass (My is relatively close to unity). For the porous 
spheres, this is slightly worse (M p is close to 2 in many cases), 
but still not very extreme. This is a consequence of the fact 
that the particles we consider are not extremely large. For very 
large particles, ry > 20 fim, the error in the mass estimate will 
be much larger. For the intermediate fractal dimensions and for 
the irregularly shaped compact particles we make an error in 
the mass estimate of 30% for the largest values of ry when us- 
ing homogeneous spheres to fit the spectra. 

In several stu dies of the 10 um spectra of circumstel- 
lar disks (see e.g. iBouwman et all 1200 it iHonda et"ail |2004; 
Ivan Boekel et all l2004 l2005t) particle size effects have been 
modeled by using a typical size for the large particles of 




2 4 6 

r, H 



Fig. 5. The best fit size of a homogeneous sphere to the 10 fim 
amorphous silicate spectra, rf, of particles with different shapes 
as a function of the volume equivalent radius of the particles, 
ry. The top panel is for the Gaussian random spheres with dif- 
ferent values of cr, while the lower panel is for the fractal ag- 
gregates with different fractal dimensions, Df. The values of <x 
and Df are indicated at the right side of the curves. The dotted 
lines represent the case in which the fitted size and the real size 
are equal, i.e. r v = rf. 



ry = 1.5 or 2.0 fim. The spectral signature of particles with 
these sizes is reported in these spectra. Since this grain size 
is significantly lar ger than that derived for the dust in the in- 
terstellar medium jKemper et all 120041) this is interpreted as 
grain growth. In these studies the particles are considered to be 
compact, and in most cases homogeneous spheres. However, 
models of grain growth in astronomical environments show 
that aggregates with a wide range of fractal dimensions can 
be formed. As we showed, in general, fractal aggregates and 
irregularly shaped compact particles display the spectral signa- 
ture of smaller homogeneous spheres. Thus, when interpreting 
emission spectra of astronomical objects using compact parti- 
cles, the size of the emitting grains is generally underestimated. 
A good estimate of the real particle shape or fractal dimension 
may be available from, for example, theoretical arguments. In 
that case the r/ derived from fitting 10/vm spectra of astronom- 
ical sources using homogeneous spheres can be used, in com- 
bination with Tabled to estimate the volume equivalent radius 
of the particles, ry. 

For example, say that we derive from the shape of an ob- 
served lOjum feature using homogeneous spheres that the typ- 
ical particle size rf = 1.7 fim. Furthermore, we expect that 
the grains present in the environment we are observing are 
fractal aggregates with Df « 1.8. We can then estimate from 
Tabled that the true volume equivalent radius of the grains is 
ry ss 6 fim. The linear extent of the aggregates can be estimated 
using y to be r c w 36 fim (see Fig.©. In addition, from Table|2] 
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Table 1. The radius of the best fit homogeneous sphere, rj (left columns), and the radius of the best fit porous sphere, r p (right 
columns), in /im for the absorption spectra of the Gaussian random spheres and the fractal aggregates. 
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Table 2. The mass scaling factor of the best fit homogeneous sphere, Mf (left columns), and of the best fit porous sphere, M p 
(right columns), for the absorption spectra of the Gaussian random spheres and the fractal aggregates. 



it can be found that the derived dust mass is a likely underesti- 
mate of the real dust mass by roughly 20%. 

Since the curves for the low fractal dimension are very flat 
for rv > 1 pm, the estimate of rv strongly depends on the ac- 
curacy of rj. Furthermore, the fractal dimension that is em- 
ployed can change the estimated value of ry from a given value 
of rf significantly. Coagulation calculations and experiments 
play a crucial role in determining the possible shapes and frac- 
tal dimensions that can be formed under different conditions. 
Therefore, these calculations and experiments are very impor- 
tant for a correct analysis of 10 pm spectra. 

Models of grain growth predict that, for example, in pro- 
topl anetary disks, grain growth rapidly removes small grains 
(Dul lemond & DominikL l2005h . When homogeneous spheres 
are used to compute the 10 pm spectra of amorphous silicate 
particles, already for a grain with a 4 pm radius the feature 
is flattened so much that it becomes very hard to detect. This 
contradicts observations of protoplanetary disk s where promi - 
nent 10 /urn features are found (see e.g. Bo uwman et al.U 2001 ). 
iDullemond & Dominikl J2005b propose a possible explanation 
for this discrepancy by considering destructive particle colli- 
sions. This causes an equilibrium situation with a significant 
amount of small grains. In addition to this, we have shown that 
when the particles grow as fluffy aggregates, the 10 pm feature 
will be visible for much larger particle sizes. 

A flattened, square 10 pm feature as seen in the ab- 
sorption spectra of homogeneous spheres is observed in the 
lO^m emission spectr a of some protoplanetary disks (see e.g 
iBouwman et all l200lh . Since the features of very irregularly 
shaped particles or fluffy aggregates as presented in Figs.[3]and 
|4]show a somewhat more rounded top, one could argue that this 



suggests that the particles in these protoplanetary disks may be 
relatively compact. However, a firmer analysis of the observa- 
tions is needed in order to confirm this. 

6. Conclusions 

We have presented calculations of the 10pm absorption spec- 
tra of complex shaped compact particles and fractal aggregates 
with various fractal dimensions. We have compared the result- 
ing spectra and studied the dependence of the spectral signature 
on the size of the particles. It is clear that the size dependence 
of the spectral signature of fractal aggregates depends strongly 
on the fractal dimension of the aggregate. For very fluffy ag- 
gregates, i.e. low fractal dimensions, the spectral signature of 
very large aggregates still looks like that of very small com- 
pact particles. We have investigated if the absorption spectra of 
complex particles can be approximated by those obtained using 
volume equivalent porous sp heres with the same volume filling 
factor as proposed bv lHage & Gre enberg ( 1990). We conclude 
that the interactions between the volume elements making up 
the particle are underestimated when we use this porous sphere 
approximation, applying Garnett effective medium theory, thus 
underestimating the effects of particle size. 

For the complex shaped compact particles the size depen- 
dence of the spectral signature also depends on the actual shape 
of the emitting grains. In general we find that the strength of the 
10 pm absorption spectrum of a nonspherical particle is equal 
to that of a smaller homogeneous spherical particle. 

Observed 10pm emission and absorption spectra are of- 
ten interpreted using compact, and in most cases homogeneous 
spherical, particles. We show that this leads to an underestimate 
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of the actual grain size when the emitting grains are either non- 
spherical compact particles or fractal aggregates. We present a 
way of estimating the true particle size and mass, when an es- 
timate of the shape or fractal dimension of the particles can be 
provided. We also show that when equivalent porous spheres 
are used to fit observed 10 //m spectra, the particle sizes are 
severely overestimated. 

Analysis of the lOjum emission spectra of circumstellar 
disks show that the spectral signature of homogeneous spher- 
ical particles with radii of approximately 2//m is present in 
these disks. When the emitting grains are actually, for instance, 
fractal aggregates with a fractal dimension of approximately 
1.8 this implies a volume equivalent radius of these aggregates 
of at least r v = 6 //m. The radius of the circumscribed sphere 
is then even of the order of r c « 36 fim. In environments where 
aggregates with even lower fractal dimensions can grow, grain 
sizes will be even more severely underestimated when using 
compact particles. 

Acknowledgements. We are grateful to the referee, N. V. 
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